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Abstract — The paper presents the effectiveness of the 
parametric integral equation system (PIES) in solving 3D 
boundary problems defined by Navier-Lame equations on the 
example of the polygonal boundary and in comparison with 
the boundary element method (BEM). Analysis was performed 
using the commercial software "BEASY" which bases on the 
BEM method, whilst in the case of PIES using an authors 
software. On the basis of two examples the number of data 
necessary to modeling of the boundary geometry, the number 
of solved algebraic equations and the stability and accuracy 
of the obtained numerical results were compared. 

Index Terms — computer modeling, computer simulation, 
boundary problems, parametric integral equation system 
(PIES), Navier-Lame equation 

I. Introduction 

The most popular in the computer simulation of three- 
dimensional boundary problems are the finite (FEM) [1,13] 
and boundary (BEM) [1,2,3] element methods. Their 
popularity comes from the variety of applications, but also 
from the wide availability of tested and therefore reliable 
software [4] . Both methods are characterized by the necessity 
of discretization at the stage of modeling the shape: in FEM 
of the boundary and the area by finite elements, in BEM 
only of the boundary using boundary elements. On the one 
hand this is advantageous, because using any number of 
such elements varied and complex shapes can be modeled, 
however, particularly in the case of 3D issues it is very time 
consuming and complicated. Another disadvantage is that 
obtained by these methods solutions are discrete, except 
that every time we get the number of solutions from boundary 
nodes (BEM, FEM) or from area (FEM), which are useless to 
us. Our previous research on solving boundary problems 
resulted in the creation and development of the parametric 
integral equation system (PIES) [7,8,9,10,11,12], which 
eliminates disadvantages mentioned above. It is an analytical 
modification of the boundary integral equation (BIE), which 
is solved using BEM. BEM has already substantially reduced 
the complexity of the boundary problems solving, because 
of the limitation of the discretization only to the boundary, so 
the next step is to completely get rid of the discretization. 
This has been achieved by analytical including the shape of 
the boundary in the mathematical formalism of BIE [7] . That 
shape, thanks to the separation of the boundary geometry 
approximation from the approximation of boundary functions, 
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can be defined in any manner. The most effective is the use 
of parametric curves in 2D problems and parametric surfaces 
in the case of three-dimensional problems [5,6]. The 
application of surfaces (e.g. Coons or Bezier) for 3D issues 
eliminates the need for the discretization of individual faces 
of the three-dimensional geometry, because they can be 
globally modeled with a single surface. The accuracy of PIES 
solutions no longer depends on the method of modeling the 
shape (if we assume that we model the exact shape using 
surface patches), only on the effectiveness of the method 
used for the approximation of boundary functions. Solutions 
in PIES are presented in the form of approximation series with 
any number of coefficients, which however means that they 
are available in the continuous version. PIES has been widely 
tested, taking into account the various problems modeled by 
2D differential equations (Laplace's, Navier-Lame, Helmholtz) 
[7, 10,9], but also modeled using three-dimensional Laplace 
[8] and Helmholtz [11] equations. The results of these tests 
were satisfactory, and therefore the study were extended on 
more complex issues, namely the problems of elasticity. 
Preliminary analysis of such problems on the example of 
polygonal issues was presented in [12]. The purpose of this 
study is to examine the effectiveness of the numerical 
implementation of PIES for solving Navier-Lame equations 
in the 3D area compared with the classical BEM. The authors 
also tried to answer the question: can PIES be a potential 
alternative to classical element method. For the implementation 
of BEM the commercial software "BEASY" was used, whilst 
in the case of PIES the authors software. Compared were: the 
number of input data needed to define the shape of the 
boundary, the number of solved algebraic equations and the 
accuracy, stability and reliability of results. Presented tests 
are beginning of research in this direction and are limited to a 
very elementary form of the area. 

II. PIES in 3D Elasticity Problems 

PIES is an analytic modification of the classical boundary 
integral equations (BIE), which are characterized by the fact 
that the boundary is defined generally by using the boundary 
integral. BEM used to solve BIE is characterized by the 
discretization of the boundary, as physical and effective way 
of its definition. Such definition of the boundary has also 
some disadvantages, because any minor change requires a 
re-discretization. Analytical modification of BIE consisted in 
defining the boundary not using the boundary integral, but 
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with surfaces used in computer graphics. These patches have 
been integrated into kernels of obtained PIES. PIES therefore 
is not defined on the boundary just as BIE, but on the 
parametric reference surface. Such a definition of PIES has 
the advantage that changing the shape of the boundary 
causes automatic adjustment of PIES to the modified shape 
and does not require classical discretization. The way of 
obtaining PIES for the Navier-Lame equations in 3D areas is 
a generalization of the technique discussed in details in [10] 
for 2D problems. PIES for considered in the paper problems 
is presented in the following form [12] 



(1) 



0.5m, (Vl h0 = j£ J J {JJ, " (v, , w>, v, w)p j (v, w) - 

P* (v, , w;, v, w)m . (v, w) }/ (v, w)dvdw, 

where v,_, < v, < v,, w,_, < w l < w n v H < v < v. , 

w^, < w < w. , / = 1,2,3..., n and « is the number of 
surfaces which define 3D geometry. 

Integrands U* (v, , w, , v, w) , P j ,*(v 1 ,H' 1 ,v ) H') from(l)are 
presented in the following matrix form 

~u„ u u 

U ,j'(y x ,w x ,v, w) = 



1 



16;r(l-v)//77 



U 2 

u. 



u 2 
u. 



Ph"{y x , w,,v, w) = 



8^(1 -v)?7 2 





Pn 


^3 




P* 


^23 




P 

1 32 


P 



(2) 



(3) 



The individual elements in the matrix (2) in the explicit form 
can be presented 



V 



Mi 
n 1 



u =^ u 



= (3-4v) + 



^2 



^23 = 



Ml 
^ 2 



C/ 33 =(3-4v) + ^ ) 
77 ^ 



whilst in the matrix (3) are presented by 

^,=((l-2v) + 34)| !Z -, J P 22 =(d-2v) + 34)?-, 
t] on t] on 

P 33 =((l-2v) + 34)^-, 
77 on 

^ = 3 ^7i^72 g'? (1 2v ^V 1 n 2 -T] 2 n, 



77" 9n 77 

7, 3 = 3 ^_ ( i_2v)Mi 
77 dn 77 
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p =3 V 2 V I (1 2v ) T MhZlh!h 
rf dn 77 ' 

p _ 3 ri 2 ?],dri 2v ^ 2 n, -rj,n 2 
rf dn 77 

p =3 '?3 ? ? 1 ^^7 (1 2v) ? ?3" 1 - ? 7 1 "3 
rf dn 77 ' 

rf dn V 
where function /^(v.w) is the Jacobian, n x ,n 2 ,n 2 are 
components of the normal vector « ( to the segment of the 
boundary marked by J , whilst 

7 7l = P« (v, w) - P/" (V, , W, ) , 77 2 = /f (v, w) - P/ 2) (V, , w, ) , 
773 =Pf(v,w)-P/ 3) (v 1 ,w 1 ) and r, = [rf +r, 2 2 +r, 2 , f , 

where P- !) (v, w), pj 2) (v, w), pj 3) (v, w) are the scalar 
components of the vector surface 
p. ( v , W ) = [P/ 11 (v, w), P/ 2) (v, w), P <3) (v, w)] T which depends 
on (v, w) parameters. This clause is also valid for the 
surface marked by / with parameters v, w, , i.e. for j = I 
and parameters v = v, and w = w 1 . 

III. Numerical modeling of 3D elasticity problems by 
PIES 

In view of the separation in PIES the approximation of the 
boundary geometry from boundary functions we deal with 
perfecting both these approximations independently. Thus, 
we use the way of the modeling of the boundary geometry - 
taken from computer graphics - which consists in modeling 
the boundary by surfaces. Such modeling requires the least 
number of the necessary data. 

A. Modeling of the three-dimensional boundary geometry 

The proposed in the paper method is characterized by the 
fact that in the mathematical formalism of PIES the boundary 
geometry is directly included, and it can be physically defined 
by a suitable connection (arbitrary) of parametric surfaces. 
This strategy and its effectiveness has been tested 
considering issues defined by the Laplace and Hemholtza 
equations [8,11]. Both polygonal and curvilinear shapes were 
tested, to define which were used respectively Coons and 
Bezier surfaces. This paper, in connection with the initial stage 
of research in solving 3D problems defined by Navier-Lame 
equations, is limited to examining polygonal areas, hence the 
description of modeling the boundary geometry will concern 
only such cases. Modeling of polygonal areas in PIES is 
performed in very simple, intuitive and effective way. Each 
face of the 3D area is defined by the Coons surface, which 
require posing only four its corner points. It is a minimal set 
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of data necessary to the accurate definition of considered 
shape. The strategy of the global modeling in PIES on the 
example of the cube with a comparison to BEM is presented 
in Fig. 1. 




a b 

Figure 1. Modeling of the 3D boundary geometry: 
a) in BEM, b) in PIES 

Presented in Fig.l cube was modeled: in BEM (Fig. la) by 
discretization of its faces using together 192 triangular 
boundary elements, whilst in PIES (Fig. lb) using only 8 corner 
points which define 6 Coons surfaces which model cube 
faces. On the basis of presented figure the efficiency of the 
proposed modeling technique is clearly seen, and manifests 
itself in reduced workload related to the modeling of the 
boundary geometry and eliminating the need for discretization. 

B. Approximation of boundary functions 

Besides the modeling of the boundary geometry in PIES 
the second problem is its solving. As it was repeatedly 
emphasized in other papers [7,8,9,10,11], the solution is 
practically reduced to the approximation of boundary 
functions on individual surfaces. For that reason we use 
approximation series with arbitrary base functions (for 2D 
problems Chebyshev, Laguerre, Legendre and Hermite 
polynomials were tested). For the each surface which model 

j segment of the boundary they are presented in the 
following form [11] 

uXv^)=YT «->r«(v)r-(w) : 

where u)""^)"" are unknown coefficients, Tj"\v), T) r) {w) 

are Chebyshev polynomials of the first kind, whilst N ,M - 
are numbers of these coefficients on the surface in two 
directions (together k = N xM )• 

Such technique is much more effective then re- 
discretization in BEM, because does not require modification 
of the defined geometry and is reduced to manipulation of 
two parameters N and m m the computer software. After 
solving PIES only solutions at the boundary are obtained. 
To receive the results also in the area the integral identity 
should be used (Appendix A). The method of receipt is the 
same as in the case of 2D problems, presented in [10]. 
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In Appendix B, are also included formulas used to calculate 
stresses in the authors software with integrands in explicit 
form. 

IV. Analysis of Results 

In the first example we consider the unit cube presented 
in Fig. 1 . One face of the cube is firmly fixed and the opposite 

is subjected to the surface compressive force p_ = IMPa . 
Material constants adopted for the calculation are: the 
Young's modulus E = XMPa and the Poisson's ratio v = 0.3 • 
In order to define the shape of the boundary in PIES we use 
rectangular Coons surfaces. Their definition is very effective 
and intuitive and is limited to posing the coordinates of the 
corner points of each surface that creates the boundary of 
the cube. As shown in Fig. lb, modeling of the boundary 
geometry leads to the use of six Coons surfaces and posing 
8 corner points of these surfaces. In the case of BEM action 
would be much more laborious, because would require the 
discretization of the boundary by boundary elements, whose 
number would depend on the desired accuracy of the 
solutions. PIES is characterized by the fact of the separation 
of the approximation of the boundary geometry from 
boundary functions and, therefore, modeling is carried out in 
a very effective way, because with minimal number of data 
required to the accurate mapping of the shape. Thus, once 
defined geometry is constant throughout the process of 
numerical calculations. The accuracy of the solutions in PIES 
depends on the number of expressions in series which 
approximate boundary functions (4a, b). It is easy to change 
parameter k , and the values used for the solution of the 
analyzed example, are as follows (for each surface separately): 
16, 25, 36, 49, 64, 81 and 100. The greater value the more 
algebraic equations we have finally to be solved, but also the 
more accurate solutions. Thus, when we choose value 100, 
to obtain the final results we must solve the system with 1800 
algebraic equations. The same example is also solved using 
the "BEASY" [4], representing the BEM. In order to analyze 
the accuracy and stability of solutions the boundary 
geometry was modeled using a different number and various 
kinds of boundary elements. Therefore, we have used 
constant, linear and quadratic boundary elements in the total 
number of 24, 54, 96, 150, 294, 600 and 864. The stabilization 
of the results (no change to third decimal place) was obtained 
using linear and square elements in the number of 600 and 
more. Assuming that we take into account 600 square 
elements, it is necessary to solve 6138 algebraic equations. It 
is over 3 times more than in the method developed by the 
authors. Table I shows the number of data needed to define 
the shape from Fig. 1 using two mentioned methods by which 
the obtained solutions are stable. The comparison of the 
number of data is more beneficial for PIES. It comes from the 
separation of the approximation of the boundary geometry 
from boundary functions and from the fact that we use the 
efficient way of the modeling of the boundary and the 
effective method for the approximation of boundary functions. 
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TABLE I. 

COMPARISON OF THE EFFECTIVENESS OF MODELING IN BEM AND PIES 





BEM 


PIES 


Numb 2i and kind of 
data for the b oundarv 
Eeometry definition 


600 quadratic 
boundary elements 
and 5400 nodes 


E comer points of 
surfac es 


Number ot solved 
alEebraic equations 


61 3B 


:soo 


Obtained results 


in3069 nodes on the 

b oundary and 
arbitrary p oints in the 
domain 


in arbitrary p oints 
on the boundary 
and in the domain 



Despite the benefits of PIES in the reduced number of data 
for modeling and the number of algebraic equations in the 
system of equations we must consider especially the 
accuracy of solutions. For this reason we analyze 

displacements u„ and stresses <j x ,<j_ in two cross-sections 

of the considered area: x = 0, y = 0,-0.5 < z < 0.5 and 

x = 0, -0.5 < y < 0.5, z = 0.3 . The results obtained by 
comparison of PIES solutions with solutions achieved using 
BEM ("BEASY") are presented in Table II and Fig. 2. 

TABLE II. 

Comparison of displacements U_ obtained by BEM and PIES 



s 


y 


z 


BEM 


PIES 


0.0 


0.0 


-0.4 


-0.07167 


-0.072275 


0.0 


0.0 


-03 


-0.15635 


-0.157760 


0.0 


0.0 


-02 


-0.250E2 


-0.253072 


0.0 


0.0 


-01 


-0.35114 


-0.354127 


0.0 


0.0 


0.0 


■0.45406 


-0.45761E 


0.0 


0.0 


1 


-0.55754 


-0.561490 


0.0 


0.0 


0.2 


-0.66053 


-0.664752 


0.0 


0.0 


0.3 


-0.76271 


-0.767112 


0.0 


0.0 


0.4 


-O.B6420 


-o. Beside 


0.0 


-0.4 


0.3 


-0.769^6 


-0.774406 


0.0 


-03 


0.3 


-0.76755 


-0.772130 


0.0 


^02 


0.3 


-0.76515 


-0.76964E 


0.0 


-01 


0.3 


-0.76336 


-0.^67793 


0.0 


0.1 


0.3 


-0.76336 


-0.767793 


0.0 


0.2 


0.3 


-0.76515 


-0.76964S 


0.0 


0.3 


0.3 


-0.76755 


-0.772130 


0.0 


0.4 


0.3 


-0.76976 


-0.774406 



As shown in Table II and Fig. 2 obtained by the proposed 
method results are very similar to the solutions obtained using 
BEM. It should be emphasized, however, that the results 
obtained by the authors require less workload related to the 
modeling (a small set of corner points) and numerical solution 
(over 3 times less number of algebraic equations in the 
equation system). Fig. 3 shows how the values of 

displacements u„ in the chosen point (0,0,-0.4) become 

stable, taking into account both methods and different values 
of parameters responsible for improving the accuracy of the 
solutions (eventually expressed as a number of algebraic 
equations). 




BEM. jj 



Figure 2. Comparison of normal stresses in cross-sections: 
a) x = 0, y = 0,-0.5 < z < 0.5 and b) x = 0,-0.5 < y < 0.5, z = 0.3 



I 1 



lllll 



inn" 



Figure 3. Displacements for different size of the algebraic equation 
system in: a) PIES and b) BEM 

As shown in Fig. 3, the stabilization of solutions with the 
desired accuracy in PIES took place at a much smaller number 
of data necessary to solve the problem than in BEM. As 
mentioned earlier it requires more than 3 times smaller system 
of algebraic equations then in the classical element method. 
In order to confirm the reliability of the proposed method and 
its effectiveness we solve the similar problem by changing 
only the boundary conditions. Thus, for the problem we 
assume the same shape and dimensions, the same material 
constants and the same degree of discretization in MEB and 
the number of expressions in approximation series in PIES. 
These assumptions caused that the data from Table 1 are 
valid also for the second considered problem. This time, 
however, the cube is subjected to the load in the form of the 

surface force p y = \MPa . 

The values of solutions obtained by BEM and PIES were 
analyzed again. One previously mentioned cross-section was 

chosen ( x = 0, y = 0,-0.5 < z < 0.5 ), and research was 
related to functions such as: displacements u y and stresses 

<J . The results are presented in Fig. 4. 

As in the previous example, obtained by PIES solutions are 
practically the same as those obtained using BEM. So we 
can consider them to be reliable, due to the fact that BEM is 
the method developed through years and widely tested. 
Simultaneously, the same level of accuracy of solutions was 
obtained with: a smaller number of data for modeling, the lack 
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of need for discretization the boundary, a smaller number of 
equations in the final system of algebraic equations and 
possibility for obtaining solutions at any point on the 
boundary and in the domain, without necessity of maintaining 
the results from all boundary nodes. 




Figure 4. Displacements and stresses obtained by BEM and PIES in 
the chosen cross-section 

Conclusions 

The paper presents the efficiency of PIES in solving three- 
dimensional polygonal boundary problems defined by Navier- 
Lame equations in comparison to the classic BEM represented 
by "BEASY" software. Obtained by PIES solutions (both 
displacements and stresses) are characterized by high 
accuracy in comparison with the element method, but also 
are obtained in a much more efficient way. It was shown that 
the proposed approach requires fewer data to model the 
boundary geometry, the system with fewer number of 
equations and unknowns is solved and does not require to 
maintain redundant solutions collected from all grid nodes, 
but it is possible to obtain them at any point of the boundary 
or area. The authors are aware of the fact that the analyzed 
example is very basic, but getting even these results require 
intensive workload. The reliability of obtained results is 
encouraging to generalize the method to more complex shapes 
of areas. 

Appendix A Displacements in a domain 

A modified integral identity for displacements takes the 
following form 



«(*) = X J 'j'(x,v,w)pj(v,w)- 

i=l iv,» v , 

Pj" {x,v,w)u j (v, w)J } (v, w)dvdw. 
Integrands from (5) are presented by 



(5) 



U "(x,v, w)- 



P"(x,v,w) = 



1 



l67T(l-v)/ur 





u a 


u 






u 


U 31 


u 32 


u 



(6) 



-1 



8;r(l-v)F 



p 

22 

R 



(7) 



Individual elements of the matrix (6) in explicit form can be 
represented 

r 2 ~ r r rr ~ r r 

t/ n =(3-4v) + ^, t/ 12 =4f, U 13 =^-,U 21 =^, 
r r r r 

"2 **** **** — ** 

w r r r ~ r r a r r 

U =(3-4v) + - 2 - U =— U =— U = — 

^ 22 W 7 '/ T -2 ' 23 -2 » 31 -2 ' 32 - 2 ) 



c/ 33 =(3-4v) + ^, 
r 

whilst from the matrix (7) take the following forms 



r dr 

^=(5 + 3^)^-. 

r on 



/>„ = (£ + 3^)— 
r on 

r 2 dr 
P 33 =(B + 3^)^, 
r On 



r 2 r ' r 2 r 



II F„ Fi„ _ 7 „ 

p _ ^ ^3^1 g r 3 W l ~ *i W 3 p _ ^ ^3^2 g r 3 n 2 ~ r 2 n 3 

r 2 r r 2 r 



dr 

where A = 3— , 5 = (l-2v) and r. = (v, w) - x , 

on 

f 2 = /f (v, w) - x 2 , f 3 = if (v, w) - x 3 , F = [jj 2 + r 2 + ~r 2 f , 
whilst x = {x 1 ,x 2 ,x 3 ] . 

Appendix B Stresses in a domain 

In order to obtain stress components we should use the 
following expression 



(T(x) = ^j j{D.'(x,v, w)p j (v, w) - 

S . * ( x , v, w) u ( v , w) } J (v, w)d vdw, 
where integrands from (8) can be expressed 



(8) 



D."[x, v, w)-- 



8w(l-v)F 2 



An 


A ,2 


A 


A*. 


A22 


A 


A 3 , 


A 3 2 


A 


An 


A,2 


b 


A 2 , 


A22 


b 


A 31 


^232 


b 


An 


A 12 


b 


A 21 


^322 


b 


A 31 


^332 


b 



(9) 
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S '{x,v, H>) = 







Sm 


L 


S m 


5,23 


s m 


^132 


^133 


"211 


^212 


^213 




^222 


^223 


^l 


^232 


^233 




^312 


^313 


L 


^322 


c 

^323 


^331 


^332 


^333 



(10) 



Elements of the matrix (9) in explicit form can be presented 



An 


= Br l + 3?/ , D m = -Br 2 + 3r 2 r 2 , 


D ii3=-- B ' 3 +3 '; 2? 3 


D m 


= Br 2 + 3r 2 r 2 , 


A 2 , = Br x + 3r x f 2 , 


A 23 = 3?;v 3 , 


K 


= Z?F 3 + 3r^r 3 , 


A 32 = 3^F 3 F 2 , Z) 133 


= Br x + 3r x r 2 , 




= Br 2 + 3r 2 r 2 , 


A, 2 =Br i +3r 2 2 r i , 


213 ^213 ' 


Dm 


= -Br 1 + 3r 2 r x 


, A22 = B h + ^ ' 


D 22i = -Br^ + 3F 2 2 F 3 


K 


— 3r r v » D 

J, 2 31 ' 232 


= BF 3 +3F, 2 F 3 , A 33 


= Br 2 + 3F.F, 2 , 


An 


= Z?F 3 + 3F 3 F, 2 , 


A12 = 3 W , A,3 


= Br x + 3r 2 r x , 


A 2 , 


= 3r r r , Z) 

J T21 ' 322 


= Z?F 3 + 3F 3 F, 2 , Z) 32 


= Br 2 + 3F 3 2 F 2 , 


K 


= -Br x + 3r 2 r x 


> A32 = + 3F, 2 


F 2 , D 333 = Br 2 + 3F 3 3 



whilst in the matrix (10) can be presented by 
S m = A(Br^ ■+ 2i?i — Sr{) + (3vnjr{ + B(3njr{ + 2jt,) — C% 
S nl = A{3F.^ — Sr-fr-) ■+ 6i^^K + S(3»^*) — Crt~ , 
i„_, = - 5Jfj|)+ 6i^r 3 + BOn^i ) - Cr^, 
S ul = -4(1^ — 5f(yQ + Sv^r-r^ — ri.r{') — Bilnf-fc +« E ), 
5]- = — 5^r 2 ") — 3i/(« 1 5i"' — raJ^f-) — 3(3«-.^f 2 — 

■Sjii = 5^F : F,)-i- 3i/(n ] ? 2 ? : i + ) — 5(3r%?}F-X 
-S B] = ^(ir 3 - 5?f/|) + 3r(«i^,?5 -1- «/, : )+ fi(3»^fj -1- 
S al = A(-5t$,t 5 ') + 3vfa?vf, + J%5{? 1 ) + J(3»,7i*|), 
- ^(wj - S^fj 1 ) + 3v(«,?i I + + 5(3ji%^H - f%J 
■Si:: = ■^('^i — 5jy^*) — 3i/(B;.^ — f^fi^) — BilnfJ^ — j%\ 
S lu = A(iTj — 5a"?5) — Si^JVl^ — Vi") — 3(3re : 'i'] — ^ X 
= 5?.^)— 3r(f%^^ — n i r z r i ') +B(5rh%r{), 

S w = A(Br 2 + 2ir ; — 5?/) -1- 6i^? ; " + ^(Sw!^ 2 + 2«,) — Crt., 
5^ = A(Br h — 5?/^) + 6mJ : 1 r l + 8(3^?^ ) - 0%, 
5-,] = A^—SP^Pj)— Iv^r^?;?} — — 5(3riji--R), 
pS'qj = — 5r ; "r,.) + 3i/(?* nfv + Hj,^") + 5(3?^^ + «,.), 
pSqj = — S^K,*)— St'fri. ^" + Ki^f,.) + 5(3Bi,FJ| + n-), 
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l 3U = J(-5K^ ) + 3F(fj : ^£ + + 5(3H,^f,X 
pSjb = -5^'f,) + 3i^H/ l ^ +H,^ ; )+ J S(3n^^ + fj ] ) : 

5 ;:] = J(-5?\f.^) + 3i- / (h ! ^ + fli^fj) + 5(3^^X 
^jjj = ^i(v^ - S^r/) + 3v(H,r : i + nJr.? : ) + BQnJs, +n,)._ 
= A(ifj - 5^)+ ivfaf/i + Hjif ) + 5(3«j^R + h,X 

S ii2 = A(Bk -5rrr : ) + 6m£r. + BQnJ^-C^ 

S w = A(Br i + 2w } -5?l) + 6wyf + BQnj? + 2«,) - Cflj, 
where C = (l-4v). 
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